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ABSTRACT 

Measuring small -signal admittance versus frtquency and forward bias 
voltage together with a new transient measurement apparently provides the most 
reliable and flexible method available for determining back surface 
recombination velocity and low-injection lifetime of the quasineutral base 
region of silicon solar cells. The new transient measurement reported here is 
called short-ci rcuit-current decay (SCCO). In this method, forward voltage 
equal to about the open-circuit or the maximum power voltage establishes 
excess holes and electrons in the junction transition region and in the 
quasineutral regions. The sudden application of a short circuit causes an 
exiting of the excess holes and electrons in the transition region within 
about ten picoseconds. From observing the slope and intercept of the 
subsequent current decay, one can determine the base lifetime and surface 
recombination velocity. The admittance measurement previously mentioned then 
enters to increase accuracy particularly for devices for which the diffusion 
length exceeds the base thickness. Detailed mathematical treatment and 
phenomenological reasoning provides support for the measurements reported and 
for the view that the method described is superior to others now existing. 

The measurment of surface recombination velocity is integrated with new 
findings about decreased values of this parameter attained via a CVD 
polysilicon/silicon interface. 
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I. INTRODUCTION 


A main purpose of this study is to establish methodologies by which one 
can experimentally determine, with good accuracy, the recombination parameters 
of silicon solar cells. These parameters include carrier lifetimes in the 
quasinuetral base region and the back surface recombination velocity including 
the velocities at internal surfaces such as those formed by a BSF region. The 
methodologies sought would fill a gap in the existing capabilities to measure 
accurately such parameters, a gap whose presence is apparently not fully 
understood by some workers in the field. The methodologies sought are to be 
flexible, in the sense that they would apply to a wide range of different 
solar-cell designs and in the sense that a subset of them would work for in- 
process control in manufacturing. Further, we seek methodologies firmly 
rooted in physical theory to avoid thereby possible misinterpretations of data 
and to provide quantitative grounds for comparisons of different experimental 
techniques. 

A review of literature will suggest the existence of some confusion 
concerning surface recombination velocities. Thus we offer remarks to aid 
understanding of the underlying physics and of the measurement of 
recombination parameters. 

Insofar as minority carriers are concerned, in a one-dimensional model 
one may regard a solar cell as a p/n junction diode bounded by front and back 
surfaces characterized by surface recombination velocities: Sf ront and 
Sb a ck* The surfaces may be free surfaces, constituted perhaps of a mixture of 
metal and thermal oxide, or internal surfaces, such as those adjoining the 
low-high junction that constitutes a back-surface field (BSF) region. This 
characterization places emphaslis on the recombination that can take place at 
the surfaces, though it includes volume recombination within the cell. 



For well over a decade, workers in photovoltaics have understood that the 
presence of a low-high junction In the BSF cell can yield an effective surface 
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recombination velocity S back ;on the low-doped side of the low-high junction 

i 

that can be orders of magnitude below the surface recombination velocity at an 
ohmic contact (which is o the order of 10® cm/s). Accurate measurement of 

i 

Sback* together with that of carrier lifetime in the quasi-enutral base, 
however, has presented problems. 

More recently the importance of passivating the front surface 
recombination velocity to increase the power conversion efficiency has become 
recognized. The first recognition of this importance, by lies, appeared in a 
final report (NASA 1974) of narrow distribution. The first full discussion of 
the importance of Sf ront in the journal literature, which emphasized 
experimental evidence in conjunction with a modeling of highly doped Si, by 
Fossum, Lindholm, and Shibib (1979), met with resistance at the outset because 
of the inertia of the dead-layer concept of Allison and Lindmayer (1973) of 
early solar-cell theory. Gradually an appreciation of the importance of 
Sf ront has emerged. One may conjecture that this emergence resulted in part 
from the understanding that the huge drift field acting on minority carriers 
in a diffused front layer that arises in customary Si p/n junction theory is 
absent because of the dependence of the Si energy gap on the shallow-level 
dopant concentration. 

Thus we recover the model stated in the first sentence of this 
INTRODUCTION: the view that a solar cell is a p/n junction bounded by front 

and back surfaces characterized by Sf ront and S back . We stress the importance 
of experimentally determining Sf ront and S back , including the determination as 
a function of the fabrication steps used in manufacturing. 
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From a theoretical viewpoint, this model yields to existing approaches 

i 

for the standard solar-cell conditions: steady state (time Independence) of 

the excitation (applied voltage or Illumination). But from an experimental 

| 

standpoint, steady-state excitation will not suffice for a widely applicable 
experimental determination of Sf ront , S^ ac | c , and other parameters (such as 
carrier lifetime) needed for Informed design. Moreover, the rapidity of 

I 

measurement by transient response of these parameters makes transient 
measurements, if accurate, attractive for In-process control at key steps In 
manufacturing. As we shall see, the use also of the real and imaginary parts 
of the small-signal admittance as a function of frequency w and forward 

voltage V provides a powerful adjunct to transient measurements. 

! 

A secondary purpose of this research is to advance the understanding of 
decreases in surface recombination velocity afforded by polysilicon/silicon 
heterojunctions at free surfaces. 

In the presentation to follow we first outline a mathematical method .hat 
systematically and compactly describes the large-signal transient and small- 
signal frequency responses of solar cells and of related devices such as 
transistors and diodes. This mathematical framework unifies the comparison of 
available methods for determining carrier recombination lifetime and surface 
recombination velocity of quasi-neutral principal regions of the devices. 

Second, exploiting this description, we survey the adequacy of various 
experimental large-signal transient methods for deducing these parameters. 
The survey is indicative, not exhaustive. 

Third, we examine in detail, both theoretically and experimentally a 
method that apparently has not been much explored previously. We demonstrate 
that this method yields both the back surface recombination velocity and the 
recombination lifetime of the' quasi-neutral base from a single transient 


I 

measurement. To Increase sensitivity and widen the range of different solar 
cells for which accurate experimental determinations can be made, we describe 
various methods that Involve measuring admittance as a function of frequency 
and forward voltage. 

Finally we briefly describe experimental research Indicating the efficacy 
of the polysilicon/sillcon heterojunction for decreasing surface recombination 
velocity. 


2. MATHEMATICAL FRAMEWORK FOR UNIFYING VIEW OF 
TRANSIENT RESPONSES AND SMALL-SIGNAL RESPONSES 

Here we develop a mathematical framework which could be applicable to 
most of the large-signal transient measurement methods and could include 
small-signal admittance methods for the determination of the lifetime and the 
back surface recombination velocity of the base region of a diode or a solar 
cell. This analysis will treat the minority-carrier density and the minority- 
carrier current in a quasi-neutral base region in low injection. Focusing on 
the quasi -neutral • base, assumed to be n-type here (of x-independent doner 
density Nq D ) with no loss in generality, will simplify the treatment. 
Extensions to the quasi-neutral emitter are straightforward, provided one 
inserts the physics relevant to m- or p* regions. 

Assume a p + /n diode in which the uniformly doped quasi -neutral base 
starts at x * 0 and has a general contact defined by arbitrary effective 
surface recombination velocity S e ^ at the far edge x = Xg N g. Such a contact 
could result, for example, from a back-surface-field (BSF) region. Assume 

l 

also low-level injection and uniform doping of the base region. Then a linear 
continuity (partial differential) equation describes the excess minority holes 
p(x,t): 
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( 1 ) 


« 

3p(x,t)/3t - O p 3 2 p(x,t)/3x 2 - p(x,t)/T p , 

where Dp Is the diffusion coefficient and Tp Is the lifetime of holes. 

If we take the Laplace transform of (1) with respect to time, we get an 
ordinary differential equation in x with parameter s: 

-p(x,0") sP(x,s) = Dp d 2 P(x,s)/dx 2 - P(x,s)/x p , (2) 

where 

P(x,s) ' / e‘ st p(x,t)dt , s ■ o+ ju, j * (-l) 1 ^ 2 (3) 

t=0 

and where p(x,0”) Is the initial condition for the excess hole density. Here, 
t = 0" denotes infitesimal negative time, and we shall treat transient 
excitation for which p(x,t) is in the steady state for t < 0. 

Thus P(x) = p(x,t), t<0 where here capital P denotes a steady-state 
excess hole density. 

Solving Eq. (2) yields 

P ( x , s ) = p(x,0”)/s + Mj exp(-x/L p ) + M 2 exp(x/L p ) (4) 

* 1/2 

where L p s (D p x p )/ ( 1 + ST p ) ' and where Mj and M 2 , given below, are to be 
determined by the boundary values at the two edges of the quasi-neutral base 
region: P(0,s) at x a 0, and p ( x qnb» s ' at xaX QNB* Substitution of (4) into (2) 
yields the steady-state continuity equation for p(x,0"), verifying that (4) Is 
the solution of (2). 

Because of quasi-neutrality and low Injection, the minority hole 
diffusion current dominates in determining the response from the quasi-neutral 
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base. The following matrix describes the density of this current at x ■ 0 and 
x a X QNB : 


’ 1(0. s) -1(0,0') 

fiO 

I 

I 


X 

1 ( X gNB * s ) '^ X QNB’ 0 ) 

*v 

e' X QNB /L p* 

e X QNB /L p* 



where i (0,0") and i( x g NtJ .0‘) are the Initial values (at t » 0") of the 
minority nole diffusion current at x-0 and x»Xg N g. In (5), nole current 
entering the quasi -neutral base is positive, by definition. 

Regarding the minority carrier densities at the two edges as the 
excitation terms for a system analogous to a linear two-port network of 
circuit theory, we have the following two-port network matrix from (4) and (5) 

m 

for the two excitations (densities) and the two responses (currents) 


i(o.s) -i(o,o‘;/s 

3 

' A n 

A 12 


_ P(0,s) -p(0,0‘)/s 

r X QNB ,S ' "'(’Vs’ 0 


A 21 

A 22 


P ( X QNB * S ) " p ( X QNB’° 


( 6 ) 


where p(0,0") and p(Xg NB .0~) are the initial values of the excess hole 
densities. Equation (6) extends a similo r earlier development [lj by 
including initial conditions so that transients may be directly studied. We 
call Eq. (6) the master equation for the quasi -neutral base, and the square A 
matrix is the characteristic matrix of the base region. In (6), ^\z sf ^z\ s 

-e( J p/ L p *) c °sech(Xg NB /L p *) and and = A 22 = e(D p /LJ). Fig. I displays the 
master equation, where the initial values are included in 1(0, s). -P(Xg NB ,s), 
etc. for compactness of expression. 

Transient solutions c u n be derived from (6) by inserting proper boundary 

% 

conditions, initial values and constraints imposed by the external circuit. 
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For example, l(0,s)*0 in OCVO open-cl rcult-voltage-decay [2], l(0,s)*constant 
for reverse step recovery [3], land P(0,s) ■ 0 for short-circuit current decay, 

the latter of which Is developed In detail here. 

| 

For small -signal methods [4]-[5], where dl, for example, Is an 
Incremental change of current, I(0,s)-I DC /s ♦ dH0,s) and P(0,s)»P DC /s ♦ 
(edV/kT)pQ C . Here the suffix DC denotes a dc steady-state variable. In later 

i 

sections, we will show briefly how to get solutions from the master equation 
for various of these methods. 

In a solar cell, the back cont' 't system, Including the low-high 
junction, is generally characterized In terms of effective recombination 
velocity , S e ff. The boundary condition at the back contact Is 
I(Xqnb* s ) “ " eS e ff P ^ X QNB* s ^* From a c ^ rcu1t viewpoint, this relation Is 
equivalent to terminating Fig. 1 by a resistor of appropriate value depenoent 
partly on S eff . Because S ef f in part determines the transient In the various 
methods named above, we can determine S e ff from the transient response, as 
will be r own. 

In the treatment to follow, we consider the utility of the master 
equation In characterizing selected measurement methods. The main emphasis 
will be placed on the short-circuit current decay. 

Before doing this, however, we shall remarl' on the simplicity provided by 
the master equation (Eq. 6) by comparing It with its counterpart in the steady 
state. 


3. TRANSIENT VS. STEADY-STATE ANALYSIS VIA TWO-PORT TECHNIQUES 
In general, the current (current density for a unit area) is the sum of 
the hole current, the electron current and the displacement current. For the 
quasi -neutral regions under study using the two-port technique described 
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previously, the displacement current Is negligible. In the steady state, the 
two-port description leading to the master equation simplifies because then 
the hole current in our example of Sec. 2 depends only on position x. This 
x-dependence results from volume recombination (relating to the minority- 
carrier lifetime) and effective surface recombination (relating to the 
effective surface recombination velocity). A two-port formulation for the 

steady state leads to the same matrix description as that derived previously, 
in which the matrix elements A^j(s) of Eq. (6) still hold but with the 
simplification that s = 0. From such a master equation, one can determine the 
hole current at the two edges of the quasi-neutral base; and, using quasi- 
neutrality together with knowledge of the steady-state currents in the 

junction space-charge region and in the p + quasi-neutral emitter region, one 
can thus find the steady-state current flowing in the external circuit or the 
voltage at the terminals of the diode. If the quasi-neutral base is the 

principal region of the device, in the sense that it contributes dominantly to 
the current or voltage at the diode terminals, then one has no need to 
consider the current components from the other two regions. 

In contrast the general time-varying mode of operation leads to a 
minority hole current in the n-type quasi-neutral base of our example that 
depends on two independent variables, x and t. The time dependence results 

i 

because the holes not only recombine within the region and at its surface, but 

also their number stored within the base varies with time. This may be 

i 

regarded as resulting from the charging or discharging hole current associated 
with 3p/3t in the hole continuity equation. This charging or discharging 

I 

current complicates the variation of the hole current in space and time. But 
the use of the Laplace transform of the two-port technique in effect reduces 

I 

the complexity of the differen 4 ”'?! equation to the level of that describing 
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the steady state; the dependence on variable t vanishes, reducing the partial 
differential equation to an ordinary differential equation in x, just as in 
the steady state. 

This comparison also brings out another point. Just as in the steady 
state, one must interpret the transient voltage and current at the diode 
terminals as resulting not only from the quasi -neutral base but also from the 
junction space-charge region and the quasi -neutral emitter. In the 
interpretation of experiments to follow, we shall discuss complications 
arising from this mul ti -regional dependence . 


4. OPEN CIRCUIT VOLTAGE DECAY (OCVD) AND REVERSE STEP RECOVERY (RSK) 

In OCVD, a widely used method [2], the free carriers in the junction 
space-charge regio.; enter to contribute to the transient. But, consistently 
with Sec. 2, and with most common usage, we concentrate at first on the n-type 
quasi -neutral base. 

From the master equation [Eq. (6)], the transient solution for the 
junction voltage is obtained from open-circuit constraint (for t > 0) that 
I (0,s )=0: 


i(0,0')L* 

p(o,s) = P (o,o~)/s 


1 ♦ {D p Ccoth(X ONa /L o *)]/L o *S eff > 
coth(X QNB /y) ♦ lD p /L p Is eff ) 


(7) 


Here we have assumed that the quasi -neutral base is the principal region in 

the sense described in Sec. 4^ that is, we neglect contributions from all the 

! 

I 

other regions of the device. 

i 

Using the Cauchy's residue theorem, we find the inverse transform of 

I 

Eq. (7): 


II 


I 


p(0,t) 


• 2i(0,0 )L p tjl + (DpX i /LpS eff )cot(X QNB K i /L p )]exp(s i t) 
1=1 eD P s 1 T P C ( i QNB / L p ) cos ec 2 ( X qnbKi ; L p ) ♦ (D p /L p S eff )] 


( 8 ) 


where is the ith singularity point (ith mode) which satisfies the 
Eigenvalue equation. 


coth(X QNB /Fv^/Lp) ♦ D p /PT7T/L p S eff = 0 (9) 

i 

i 

and K. = /-1-s.x > 0, where s. < 0. 

l l p l 

As can be seen in Eq. (8), the decay of the excess hole density at x=0 is 
a sum of exponentials; each Eigenvalue s^ is called a mode, as in the 
electromagnetic theory. Appendix A treats the details of determining the 
Eigenvalues s^ from Eq. (9) (and from the similar Eq. (11) derived below). 

The decaying time constant -1/sj, of the first mode is much the largest 
of the modes. Both Sj and the initial amplitude of the first mode are 
functions of S e ff and tp. Thus separating the first mode from the observed 
junction voltage decay curve, by identifying the linear portion of v(t), will 
enable, in principle, determination of S e ff and Tp simultaneously. But our 
recent experience, . coupled with that cited in [6], suggests that this is 
seldom possible in practice for Si devices at T * 300 K. In Si devices the 
open-voltage decay curve is usually bent up or bent down because of 
discharging and recombination within the space-charge region. 

As mentioned in [4], the mobile charge within the space-charge region 
contributes significantly to the observed voltage transient for Si, in which 
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~ 10 cm , but not in Ge, for which OCVD was first developed, and for 
13 -3 

which n.. ~ 10 cm , Here n^ is the Intrinsic density and is also the ratio 
of the pre-exponential factors that govern contributions from the quasi - 
neutral regions relative to those from the junction space-charge region. 

Thus we identify the transient decay of mobile electrons and holes within 
the p/n junction space-charge region, which persists throughout the open- 
circuit voltage decay (OCVD), as a mechanism that distors OCVD so 
significantly that the conventional treatment of OCVD will not reliably 
determine Xp or S E pp. The conventional treatment is consistent with that 
proceeding from the master equation, as described in this section. The 
interested reader may consult Ref. 4 for experimental comparisons that lead to 
this conclusion. We shall not pause here to present these. 

Rather we shall turn briefly to possible methods to remove the effects of 
this distortion. In an attempt to characterize the space-charge-region 

contribution to the observed transient voltage [6], quasi-static 

approximations and a description of the forward-voltage capacitance of the 
space-charge region based on the depletion approximation were combined to give 
rough estimates of this contribution. We plan to refine the approximations 
and the estimates in a future publication, leading possibly to a variant of 

t 

OCVD useful for determining tp and S e ff. 

In the reverse step recovery (RSR) method [3], in which again the diode 
is subjected to steady forward voltage for t > 0, we have two constraints (for 

I 

t > 0). The first is 1(0, s) = constant (reverse current), 0 < t < x s , where 

x. is the time needed for the excess hole density p(0,t) to vanish. This is 

! 2 

the primary constraint. (The second constraint is p(0,t) * for 

x s < t < ®, a result of applied reverse bias through a resistor.) 
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The primarily observablje storage time x $ Is estimated by following a 

I 

procedure similar to that djescrlbed in Sec. 4, proceeding from the master 
equation. 

This method suffers difficulties similar to that of the OCVD method. 
Because p(0,t) > 0 for 0 < t < x s , the decay of mobile hole and electron 
concentrations in the p/n junction space-charge region complicates the 
interpretation of the measured x in terms of the desired parameters, x and 

J U 

s eff * 

In addition to this, during the portion of recovery transient occuring 
for x,. < t < ®, the reverse generation current is often large enough to 
saturate the recovery current so quickly that we have no sizable linear 

I 

portion of the first-mode curve on a plot of ln[i(t)] vs. t. This linear 
portion provides interpretable data for Ge devices [3], but not often for Si 

devices according to our experiments. 

» 

We shall not pause here to present experimental evidence, but rather 
postpone that presentation until a planned future publication where comparison 

l 

with previous worl< can be comprehensive. 


5. SHORT CIRCUITED CURRENT DECAY (SCCD), A NEW METHOD: PART 1 

5.1 Brief Physics and Mathematics 

In this method, one first applies a forward bias to set up a steady-state 
condition and then suddenly applies zero bias through a resistance so small 
that the constraint is essentially that of a short circuit. Thus, for t > 0, 
the p/n junction space-charge and quasi-neutral regions discharge. One 
measures the transient current via the voltage across the small resistor. If 
the discharging time constants related to the charge stored within the quasi- 
neutral emitter and the junction space-charge region are much smaller than 
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from the quasi-neutral base,' or, e can separate the first mode of the quasi- 
neutral -base response and determine S e ff and Tp. 

We first consider the time of response of the junction space-charge 
region. Upon the removal of the forward voltage, the constraint at the 
terminals becomes essentially that of a short circuit. The majority-carrier 
quasi-Fermi levels at the two ohmic contacts immediately become coincident, 
and the junction barrier voltage rises to its height at equilibrium within the 
order of the dielectric relaxation time of the quasi-neutral regions, times 
that are of the order of no greater than 10 _1 ^s. This occurs because the 
negative change in the applied forward voltage introduces a deficit of 
majority holes near the ohmic contact of the p + emitter and a deficit of 
majority electrons near the ohmic contact in the quasi -neutral base. The 
resulting Coulomb forces cause majority carriers to rush from the edges of the 
junction barrier regions, thus causing the nearly sudden rise of the barrier 
height to its equilibrium value. (The physics governing this phenomena comes 
from Maxwell's Curl H_ = 1_ + 3D/ 3t ; taking the divergence of both sides yields 
0 = div j_ + ' 3 (di v D)/dt, which, when combined with i_ = (oD/e) and 
div£ = p, yields a response of the order of e/o, the dielectric relaxation 
time.) 

Following this readjustment of the barrier height, the excess holes and 
electrons exit the junction space-charge region within a transit time of this 
region (about 10"^s typically), where they become majority carriers in the 

. i 

t 

quasi-neutral region and thus exit the device within the order of a dielectric 
relaxation time. j 

Thus the discharging of excess holes and electrons within the junction 
space-charge region in the SQCD method occurs within a time of the order of 
lCT^s, which is much less th^n any of the times associated with discharge of 
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the quasi -neutral regions, i Thi s absence in effect of excess holes and 
electrons within the junction space-charge region greatly simplifies the 
interpretation of the observed transient. It Is one of the main advantages of 
this method of measurement . 

A more detailed discussion of the vanishing of excess holes and electrons 
within the junction space-charge region appears in Appendix B. 

The discharge of the quasi -neutral emitter depends on the energy-gap 
narrowing, the minority carrier mobility and diffusivity, the minority-carrier 
lifetime, and the effective surface recombination velocity of this region. 
For many solar cells, this discharge time will be much faster than that of the 
quasi-neutral base, and we shall assume this is so in the discussion to 
fol low. 

Having established that the mobile carriers in the junction space-charge 
region enter the short-circuit-decay transient during an interval of time too 
short to be observed, and noting also now that negligible generation or 
recombination of electrons or holes within this region will occur during the 
transient, we now turn to the observable transient current. Inserting the 
constraint, P(0,s) = 0, into the master equation, Eq. 6, leads to 


eD p(0,0~) 

I(0.s) = 1(0,0 )/s 

P 


c°th(x QNB /L D ») ♦ p D /ys eff 

1 + *- D p /L p* S eff- Jcoth ^ X QNB /L p*^ 


( 10 ) 


Cauchy's residue theorem yields the inverse transform of (10): 


cot (Kj x q NB /L p ) ~ D p K i /L p S ef l 


» eD p(0,0“)K. 

i(t) = -J — Vi 2 2 

1=1 1 p ( T p/ 2I V + ( X QNB^ 2S eff ) cosec ^ K i X QNB^ L p ^ 


s i t 

e 


( 11 ) 


i 


16 


I 



Equations (13) and (14) contain four unknowns: if-j rs t mode^* s l» T p and 
S e ^r^ . Tne parameters, Sj and i^ r$t mode (0) are determined from the 

straight-line portion of the observed decay (in Fig. 3c to be discussed 

- 2 

below): p(0,0 ) = (n 1 -/Npp)[exp(ev(0 )/kT)-l], Here v(0“) is known from the 

steady forward voltage applied for t < 0 and the doping concentration Nqq of 
the base is measured by usual methods; Dp(N DD ) is known from the standard 
tables and Xq NB is measured. Combining (14) and (13) then yields the desired 
parameters: x p and S e ff. 



I 
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Voltage across BSF #1 solar cell (vertical :.2V/div). 

Current through BSF #1 solar cell (vertical :lmA/div) . 

Log scale representation of (b) (vertical I: .lV/div) , where v(t) 

= (mkT/e)log e [i (t)/I 0 + 1]. 






We measured the voltage across the solar cell under study. As 
Illustrated in Fig. 3(a), the voltage drops by 0.1 V within 1 us. This means 
that the excess minority carrier density p(0,t) drops to 2% of Its initial 
value within 1 us. The speed is circuit limited. One could design a much 
faster circuit. Here x d = -1/Sj is the first-mode decay time, influenced by 
both volume and surface recombination in the base. But the circuit used 
suffices because ip » 1 us for the solar cells studied. Fig. 3(b) shows the 
current during the transient. Fig. 3(c) is its semi -logarithmic counterpart, 
illustrating the straight-line position of the transient obtained from the 
output of the logarithmic amplifier in Fig. 2. From this x d is determined. 
Since the voltage at node B is purely exponential for a time, the 
corresponding output voltage at node C is linear in time, as Fig. 3(c) 
illustrates. We used switching diodes in the log amplifier whose I-V 
characteristic is V = .0385 ln(I/I 0 +l). If the first-mode current is 

If i rst-mode^) = constant exp(-t/x d ) , x d = -1/Sj , (15) 

then the slope of the output voltage of log amplifier is -38.5 mV/x d . 
Extrapolation of the straight portion in Fig. 3(c) yields the initial value of 

'first-mode t 0 *) as the 

We measured the decay time constant and the initial amplitude of the 
first-mode current as follows. For DEVICE 1, x d = -1 /s ^ = 29.3 usee, i^ rst _ 
mode(0 + ) = 2.73 mA for v(0") = 0.44 V and T - 303.1 K. For DEVICE 2, x d = 

24.5 usee, ifi rst -mode( 0+ ) = 4 * 35 mA for v(0_) = °* 5 v and T 3 302.9 K - For 
DEVICE 3, x d = 28.5 usee, Virst-mode^ 0+ ^ 3 * 696 mA at v (°") 3 * 47 v and 

303.5 K. Here v(0") denotes the steady forward voltage applied across the 

I 

i 

solar cell before transient. 
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From the above development, these results give: For DEVICE 1, Tp 3 

119 ws, S e ff 3 25 cm/s; for DEVICE 2 3 tp 3 119 us, S e ff 3 60 cm/s; for 
DEVICE 3, T p 3 213 us, S e ff 3 100 cm/s. These results agree favorably with 
those obtained for the same devices by using the more time-consuming methods 
detailed in [4] - [5], 


5.3 Comparison with Other Methods 

Mo«* measurement methods for the determination of the minority-carrier 
lifetime and the surface recombination velocity of the base region of Si solar 
cells share a common problem caused by the existence of the sizable number of 
the mobile carriers within the space-charge-region. These methods, among 
open-circuit voltage decay and reverse step recovery were originally developed 
for Ge devices. Silicon has a much larger energy gap Eq than does Ge. Thus 
the distortion of the measured response by carriers stored in the junction 
space-charge region is much more pronounced in Si, mathematically because of 
the role of the intrinsic density n^ discussed in Sec. 6. 

If the electronic switch providing the short circuit closes fast enough, 
the mobile holes and electrons stored for negative time in the junction space- 
charge region play no role in determining the response of the short-circuit- 
current decay for the solar cells described earlier. In our experiments, the 
simple circuit of Fig. 2 had speed limitations, but these limitations did not 
markedly influence the accuracy of the determined base lifetime and surface 
recombination velocity fo r these solar cells (see Sec. 6 for devices for which 
this circuit must be improved). This lack of influence results because the 
decay time of the first-mode response, which accounts for vanishing of 
minority holes both by volume recombination within the quasi-neutral base and 

I 

effectively by surface recombination, greatly exceeded the time required for 
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the excess hole density at , the base edge of the space-charge region to 
decrease by two orders of magnitude. 

Apart from this potential circuit limitation, which one can overcome by 
Improved circuit design (see Sec. 6), a more basic consideration can limit the 
accuracy of the short-circuit-current decay (SCCO) method. In general, the 

current response derives from vanishing of minority carriers not only In the 
quasi-neutral base but also in the quasi-neutral emitter. For the solar cells 
explored in this study, the emitter contributes negligibly to the observed 
response because of the low doping concentration of the base and because of 
the low-injection conditions for which the response was measured. But for 
other solar cells or for higher levels of excitation, the recombination 

current of the quasi-neutral emitter can contribute significantly. 

The contribution from the quasi-neutral emitter may be viewed as an 
opportunity rather than as a limitation. That i Jt the SCCD method may have 
utility wider than that treated here. If the quasi -neutral emitter contains 
excess charge whose decay time dominates in determining the transient 
observed, in some devices the parameters of the physical electronics of the 

highly doped and thin quasi-neutral emitter can be explored using SCCD. An 

example may be a transistor with a thin and highly doped base region or a 
solar cell having a high open-circuit voltage controlled by the recombination 
current in the emitter. For such devices, the absence of contributions from 
carriers in the junction space-charge region becomes a particularly key 

advantage not offered by either open-circuit voltage decay or step reverse 
recovery. The SCCD method also may enable parameter determination if high 
injection in the base prevails. Exploring these possible uses will require 
fast switching circuits and determination of the existence of a dominant 

i 

I 

relaxation time from minority carriers in the highly doped emitter. 

I 
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Note that the SCCD method determines the base lifetime and the effective 
surface recombination velocity of a i>SF solar cell by a single transient 
measurement. One can easily automate the determination of these parameters 
from parameters directly measured from the transient by a computer program, 
and the measurement Itself may be automated. This suggests that SCCD may be 
useful for In-process control in sclar-cell manufacturing. 

Recall the mathematical formulation of the relevant boundary- value 
problem that led to a description similar to that of two-port network 
theory. The advantages of this formulation were only touched upon earlier and 
only the bare elements of its relation to open-circuit voltage decay and step 
reverse recovery were developed. Further exploitation to enable systematic 
development and comparison of small -signal and transient methods for the 
determination of material parameters of solar cells and other junction devices 
is treated. In part, In Sec. 6. 

6. IMPROVEMENT OF SHORT- CIRCUIT- CUR RE NT DECAY METHOD 
FOR DETERMINING x AND S OF THE QUASI-NEUTRAL BASE, AN IMPROVEMENT: 

PART 2 

In Sec. 5 we have proposed and illustrated a new method for measuring the 
surface recombination velocity S and the recombination lifetime x in the 
quasineutral base of silicon solar cells and solar cells of other 
semiconductors. 

In the short-circuit decay approach, for t < 0, a forward voltage has 
distributed excess holes and electrons through the volume of the solar cell. 
Then a switch establishes a short circuit across the solar cell at t - 0. 
For t > 0, the resulting transient decay Is a key observable. The decay 
results from the recombination of the excess carriers within the volume of the 

I 

device and at Its surfaces. 
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The short-circuit decay method has an advantage compared with other 
methods. This follows f r om the dielectric relaxation and the drifting of 

excess holes that yields a junction space-charge region (junction transition 
region) in which practically no excess holes and electrons exist after about 
10 ps. have passed after the switch closes. Thus this carrier storage does 
not distort the transient decay. Such distortion obscures the Interpretation 
of other transient methods for determining S and x [7]. 

For germanium, for which most of the commonly used transient methods was 
developed, the error thus introduced is negligible; the excess carrier storage 
In the junction transition region for forward voltage is much less than the 
storage in the quasineutral base region. This is not so for Si or GaAs. The 
larger energy gaps and smaller intrinsic densities of these materials lead to 
relatively very large carrier storage in the junction transition region, as 
compared with that present in Ge devices. Thus the error introduced by this 
storage for Si solar cells, particularly highly efficient solar cells in which 
the carrier storage and recombination in the quasi neutral base Is small, 
invalidates the conventional transient methods. 

6 . 1 Improvements in the Circuit for the Short-Circuit 

Current Decay 

We have improved the circuit relative to that reported in Sec. 5. An 
MUSFET switch (Fig. 4) replaces our former more complicated and slower 
switching circuit (Fig. 1 of Sec. b). Part (b) of Fig. 4 shows the 
observables for t > 0. These consist of a slope and an intercept. The 
relevant theory appears in [7]. 

The circuit of Fig. 4 offers significant advantages. Our previous 

I 

circuit led to a relatively slow establishment of the short-circuit 

I 

conditions. For certain solar cells, this made impossible the determination 
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Fig. 4(c) Experimental current decay of a n + /p/p + BSF solar cell (p^ ase = 
0.3 ftn, W hase = 367 ym, = 6.4 usee, L n = 180 ym, S e ^ ^ = 
1.3 x 10 J cm/sec). The vertical scale is 100 mA/division. 
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* of S and x. The circuit of Fig. 4 switches in approximately 10 ns., and no 

i 

problems exist in interpreting the lifetime and S of any of the over ten solar 
cells of different types we have treated. • 

6.2 Remarks on Sensitivity 

Figures 5(a) and (b), developed from theory, illustrate aspects of the 
sensitivity of the SCCD method. Let L = minority-carrier diffusion length. 
Let W = quasineutral base thickness. Then for thin cells (W << L) the method 
is sensitive to S because most minority carriers recombine at the back 
surface. Conversely, for thick cells, W >> L, the method is sensitive to x 
but not to S. 

In the limiting cases just described, the SCCD method will yield only one 
of the desired parameters, either x or S. To determine the other parameter 
requires another measurement that depends both on S and x (or L). For 
example, one can use the dark reverse saturation current or the open-circuit 
voltage. Both of these parameters may be influenced by quasineutral emitter 
recombination. Thus error can occur. Use of small-signal admittance avoids 
this error. 


6.3 Small-Signal Admittance, Preliminary Remarks 
This technique works well as the needed supplement. Neugroschel (1984) 
has described details of the application of this technique as a supplement to 
SCCD. This work resulted as part of an ongoing collaboration between the 
present Principle Investigator and Professor Neugroschel. Graduate Research 
Student, T. W. Jung, is contirtuing these efforts. We describe the method and 
its results in Sec. 7. 

I 
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Table 1: Summary of results for some typical solar cells. 

The values for L^ aS p and S e ff were obtained using 
the SSCO method, unless marked otherwise 


CELL 

p base 

(flcm) 

w base 

(urn) 

L base 

(ym) 

s eff 

(cm/sec) 


n + /p/p + 

10 

227 

454 

105 


BSF 



450* 



n + /p/p + 

10 

103 

250 

2.9xl0 3 


BSF 






n + /p/p + 

10 

360 

512 

2xl0 5 


n + /p/p + 

0.15 

295 

100 

— 


p + /n/n + 

10 

320 

oo 

°+ 

503* 


BSF 




500 f 


n + /p/p + 

10 

92 

-600* 

-180 


BSF 







HF 

* obtained from Gq N 


t obtained from G™ 
QN 
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6.4 Representative Results 

These appear in Table 1. Although over ten different solar cells were 

i 

measured using SCCD, we report six here, relegating further reporting to the 
future when the detailed physical make-up of the solar cells becomes 
avai 1 able. 


7. SMALL-SIGNAL ADMITTANCE MEASUREMENTS [9-11] 

Small -signal admittance measurements can be used to analyze a variety of 
semiconductor devices. We discuss here specifically the applications for 
analyzing the solar cells, namely measurement of the base L and S e ff and the 
separation of the emitter and the base current components. The small -signal 
measurements can be performed either at low-frequencies (wx « 1) or hign 
frequencies (wx » I). The choice of a particular frequency ranye will depend 
on the W/L ratio. 


7.1 Low-frequency method (LF) [9 , 10 ] 

Consider a n + /p/p + BSF solar cell shown ir. Fig. 6(a). For a low- 
frequency signal with ujx n « 1, where x n is the minority-carrie r electron 
lifetime in the p-type base, we derive the expressions for the small-signal 
quasi -neutral base capacitance C^g and conductance Gg^g, respectively 
(see equations (Dl) and (D2) in Appendix D) . Equations (Dl) and (D2) contain 
four unknowns: Cgjjg, Ggj^, L n , and S eff . The parameters Cg^g and Gg^g 

are measured and the combination of (Dl) and (D2) yields L n and S g ff. 

It is worthwhile to discuss in more detail a few special cases: 

7.1.1 Long diode: W > L 

For this case, (Dl) and (D2) yield a simple expression for CgJ^g and x n 
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Fig. 6 



(a) 


(a) Schematic diagram of an n + /p/p + BSF solar cell. 


* = *e + *scn + / N (x ) dx + l B « 



(b) Qualitative sketches of minority-carrier distributions in the 
dark. i 



( 16 ) 



and 


2C 


LF 


3 


NB 


(17) 


QNB 


The base diffusion length is obtained either from (2) or from (3). Thus 
2 

s ^ n T n and T n are determined independently. The details concerning the 

deduction of Cq^ b and G^ B from the data are discussed in [9]. As an 

illustrative example, we show in Tig. 7 the measured Cgj^ B (V) and Gq^ b (V) 

plots for the p + /n device with n DD = 1 .25 x 10 15 cm" 3 . The analysis using (16) 

2 

or (17) gives L p = 80 ym. Note that n i does not enter Eq. (17). Thus the 
determination of is independent of energy-gap narrowing. 


7.1.2 BSF solar cell: i L , WD /L|; < S ~ < (D„ WJ 

p n p n n ett n p 

In this case, (Dl) and (D2) are solved to yield L n and S e ff. Figure 8 

shows the measured C(V) and G(V) dependencies for a p + /n/n + BSF solar cell 

from which we derive L = 500 urn and S „ a 80 cm/sec. 

p ef f 

2 

The method fails, however, for S ~ << D W/L ; in this case 

eft p n 

GqNB * KW p^ T n ^ 1e ^ ds T n» d ut 3 eff cannot found. Another limitation exists 
for high values of S e ff >> D n /W p : in this case both (Dl) and (D2) are 
independent of L n and S e ff. 

The above difficulties with the LF method can be largely eliminated by 
the high-frequency approach. 


7.2 High-frequency method (HF) [11] 

I 

We treat the high frequency method for two special cases. 
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7.2.1 cut p > 10 and 0.1 < W p /L n < 1 

i 

The small -signal admittance then Is 


KD wt 1/2 


T- 1/2 

J"G£) ] 


The important conclusion from (18) is that the dependence gives the range 

of 0.1 i W p' ,L n ~ 1 regardless of the value of S eff • To obtain the desirable 

HF 

parameters, w? measure Gg^g vs w for ojt > 10 and extrapolate to lower 
frequencies to obtain an intercept wj with GgJ^g given by (D2). This gives 



2 

Note again that n.. does not enter. Equation (19) cannot be be solved for L n 
and S e ff except for the following cas?s: 



( 20 ) 


( 21 ) 



The method is illustrated in Fig. 9 for the p + /n/n + solar cell of Fig. 8. 
follows the .1/2 dependence for f > 10 4 Hz with the intercept at oij 3 
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u. z 



up 

Measu r ed high frequency conductance G 0( , vs freqency for the 
p + /n/n + solar cell of Fig. 8. The corrductance was measured at 
forward bias V = 0.5 V and shows nr'* dependence. 


J 



U*)10 3 l/sec. Using (19) and combining with S e ff * 80 cm/sec obtained by the 

i 

LF method, we have Lp * 503 um, which Is In excellent agreement with 

L ■ 500 um using the LF method alone. 

P 

Even though the general solution (19) cannot give L and S e ff exactly, and 
only one of the parameters Is obtained either from (20) or from (21), the 
method Is very useful because: (1) the vs dependence shows that 

L > W; (11) the vs w 1 ^ 2 dependence indicates that the emitter 

contribution to the conductance (and dark current) Is negligible (this point 
Is discussed further below); (111) the SSCD for L > W yields an accurate value 
of S e ff only and using this value in (19) we obtain an accurate value for L. 
The combination of these two methods gives L and S e ff for practically any 

cell . 

HF 

The sensitivity of the HF method to the emitter component Gq^ of the 

HF HF HF 

total measured quasi -neutral conductance Gq N * ^qnb + ^QNE is explored in 

HF 

Fig. 1C. The time constant of Gq NE is given by either the Auger lifetime 

in the heavily doped emitter, or by the combination of t a and the 

transit time [12]. The emitter time constant is much shorter than the base 

HF 

lifetime, thus Gq NE is frequency independent up to f * 1 /t^ >> 1/tg. 

Figure 10 shows G^ = G^g + Gq NE for an arbitrary choice of g qne/Gqnb« The 

region far away from the knee can be fitted to a straight line with 
G a where m > 2. Notice, however, that wj for GqJ^ >0 Is close to 

the Intercept value wj for GqJ^ = 0. Furthermore, since L a wj^, a small 
error in uj gives only a negligible error in L. For example, for 
Gq NE * 6 q N b. <*>1 “ 1.5 «! ( g q N £ * 0). This gives an error In L of only about 

ion. 
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Fig. 10 


Plot of - G^ b + gJL for different ratios of 
0, 0.5, 1,3. 



7.2.2 u)T n > 10, W p /L n i 0.1 

i 

The condition W p /L n i 0.1 may apply for the thin cells (50 - 100 urn) 
with a very long lifetime. For this case we have 




+ WD„ + 


^n + S eff W p) 


7 


S eff W p) 


( 22 ) 


For BSF cells, S e ff < D n /W p and (22) yields 

GgJ * K[» 2 (W 3 /3D n ) + S eff ] . (23) 

HF 

Figure 11 shows the Gq N vs w dependence for a 8 ym thick epitaxial 
n-type layer with doping density Nq D = 5 x 10^ cm"^. The o> 2 dependence for 
f > 1.5 MHz immediately gives L p > 10W > 80 ym and also S e ff << (u^W^/3D n << 
1.2 x 10-* cm/sec. More accurate analysis of the knee region below the uj 2 
dependence gives = 120 cm/sec and using this value in (Bl) gives more 

accurate L p = 90 ym. 

Note, that the HF method for L > 10 W gives only the lower limit of L and 
the upper limit of S e ff. A combination of this technique with either the LF 
method or the SCCD can give more accurate results. 

| 

7.2.3 u>x n ~ 10 , W p /L n ~ 0.1 

HF 1/2 

For the previous two special cases we have obtained G^g a u> ' for 

W p /L n < 0.1 and Ggj^g r. u for W p /L n < 0.1. Obviously, there has to be an 
intermediate range for W p /L n ~ 0.1 where G^g a J 11 ( < m < 2). One possible 
approach here is to obtain from the SCCD method and then fit the 

theoretical ^(S^, L n ) with the experiment. A very reasonable 
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Frequency dependence of high-frequency conductance for a thin 
(8 urn) n-±ype epitaxial layer (from Ref. [5]). The conductance 
follows or dependence. 




approximation of x n can be' made, however, by realizing that the Gg^j vs u> 

| 

dependence begins to increase from its low-frequency value for 

a) = wj = 10/-c n til]* Thus, ! x n * M/u)j, where u>j can be approximated as 
the intercept of the line with the extrapolated Gj^ a u dependence. 


7.2.4 Interpretation by Analocy. 

I 

The mathematical descriptions given above indicate that the measurements 
of Y(jw) can yield lifetime i even though the thickness W « But 

these descriptions do not fully explain this result from the viewpoint of 
physical reasoning. 

Analogy can help provide understanding of the relevant physics. Consider 
an RC transmission line having resistance r and capacitance c per unit 
length. Suppose the transmission line has no dissipative loss in parallel 
with the distributed (resistive) capacitance. Then, if v(x,t) denotes the 
electric potential and i(x,t) the current on the line, the describing partial 
differential equations are: 


a v(x,t) _ 5v(x,t) 

ax 2 9t 


5 2 i(x,t) 

r 


ax 


ai(x,t) 

at 


(24) 


Using the connections, v(x,t) -► n(s,t) and re 1/D n , we see that Eq. (24) is 
anal ogous to the continuity equation for the p-type quasineutral base region 
of a solar cell in low injection provided the lifetime i ®. Via this 
analogy, we shall now discuss the dependence on u> of the effective thickness 
of the transmission line (or the quasineutral base). 

We subject Eq. (24) to Laplace transformation with respect to t (v,(x,t) 
-*• V(x,s), i(x,t) = I(x,s)) and solve the resulting ordinary differential 
equations. This yields, for example, 
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I ( x , s ) = A exp[x/srcj + B exp[-x/srcj 


(25) 


where s is the complex frequency variable: s = o + jo>. The integration 
constants A and B are determined by the Doundary conditions, which involve 

V(x,s), and a two-port formulation results, similar to that described in 

1 

Sec. 2. This two-port description provide a convenient vehicle for 
quantitative description. 

Now consider, for example, the re transni sion line connected to a 
voltage source V^(s) at the left port and disconnected (open-circuited) at the 
right port where the voltage is V 2 ^s) Consider the transfer ratio 
V 2 (s)/Vi( s ) for the case s -► jw. Qualitatively we see that the transmission 
line is a low-pass filter. Thus az u> approaches 0 we anticipate that the 
transfer ratio will approach unity, and for large w, we anticipate that the 
transfer ratio will approach zero. To see the meaning of large o>, we 

I 

calculate the frequency for which the magnitude of the transfer ratio equals 
U.707; this is the so-called cutoff frequency for which the ratio is 3 dbs. 
below its low-frequency value. (By definition the number of decibels = 10 
log|V 2 /V 1 | 2 -). 

The two-port formulation mentioned above quickly yields the required 
condition : 

|sechL/jw o rc | = 0.707 (26) 

where ojq is the cutoff frequency. From this it follows that 
u Q 


= 2 . 43/ [ ( Lr ) ( Lc ) ] 


(27) 


in which the denominator is i the product of the total resistance and the total 
capacitance of the transmission line. Use of the analogy between the 
transmission line and the quasineutral base region of thickness W yields 

w 0 = 2.43 D n /W 2 ■ (28) 

which is the reciprocal of the transit time across the quasineutral base of 
thickness W when t -*• ®. (The more commonly seen expression in which 2 
replaces 2.43 is the quasi -static approximation for the transit time). 

Thus when the signal frequency applied to the rc transmission line, or to 
the analogous quasi -neutral base region, equals the reciprocal of the transit 
time, then the signal reaching the output port of the transmission line, or 
the back surface of the quasineutral base, falls 3 dbs below its low-frequency 
value. For the quasineutral base, this implies that the interplay of minority 

I 

electron diffusion and storage distributed through this region confines the 
input signal partially within the base region, even for the limiting condition 
that x * ® and S = 0. Hence we see qualitatively that application of a signal 
of the form exp(jwt) enables determination of x, when x is finite and large in 
the sense that the time-invariant condition W « /DT holds. This is what we 
set out to demonstrate. 

To this point we have treated (24) as an anology. In fact, it is an 
exact description for low injection conditions in the quasi -neutral base if we 
regard v(x,t) to be the electrochemical potential (quasi-Fermi potential) for 
the minority electrons and if we identify, as before, that r c -► 1/D n . That 
is, the rc transmission line exactly describes the quasineutral base for the 
conditions stipulated and for x = ®. For finite x one needs to add a 
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conductance in parallel with the distributed capacitance to attain a full 

I 

network representation of thje minority-electron transmission line. Sah bases 

I 

his numerical solutions for, solar cells and other semiconductor devices on 
generalizations of such transmission lines, the most complete versions of 
which appear in a 1971 paper [8]. 

I 

7.3 Regional Analysis of Solar Cells 
It is important to analyze the contributions of each region of the cell 
to the total dark current (or V QC ). Such an analysis is demonstrated here for 
a n + /p/p + BSF solar cell shown in Fig. 6(b). The analysis is based on the 
determination of the base parameters x n and by one of the methods 

I 

discussed earlier. This is sufficient to calculate the profile of the 
minority electrons in the base. The recombination losses in the base are 
given by (D2) and the recombination losses in the p + -BSF portion of the base 
are 


W * K * A < S eff N (V 


(29) 


The space-charge region (SCR) recombination current 1*^ can be determined 
graphically [12] and the emitter contribution is obtained by realizing that 
the total dark current is 


*D = *E + ! SCR + *8 + *B ’ 

For example, such an analysis of the p + /n/n + BSF cell of Fig. 5 gave 
[10]: L p « 500 pm, S^ * 80 cm/sec, I g » 0.8 Ip, I SCR = 0.2 Ip, I £ « I D , 
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8. COMBINED SCCD AND SMALL SIGNAL ADMITTANCE 

Table 1 repeats the summary of results for a number of different cells. 
A comparison of results obtained by different methods, shown for some cells, 
demonstrates very good agreement. Notice, in particular, the last cell in 
Table I, which is a thin cell (W base - 92 wm) with L n » W B . For this cell, 
the SCCD method gives S gff “ 180 cm/sec, but the method is Insensitive to L n 
(see Fig. 5(a)). We have to combine the SCCD method with the high-frequency 
small -signal admittance method and then use (19) with S e ff obtained from the 
SCCD method to determine L p . 

The main conclusion of this study is that the SCCD method and the small - 
signal admittance methods yield a rapid and reliable determination of the base 
parameters. They also allow the determination of the relative importance of 
the base and the emitter regions with regard to cell efficiency. 
Identification of the region limiting the efficiency is a key to an informed 
cell design. 

9. POLYSILICON/SILICON HETEROJUNCTIONS APPLIED TO SOLAR CELLS 

Replacement in Si bipolar junction transitor of the metal emitter contact 
by a highly doped polysilicon layer improves the common-emitter current gain 
of bipolar transistors. The principal investigator has reviewed some twenty 
papers dealing with this issue, prominent among which are [14], [15], [16], 
and [17]. According to de Graaff and de Groot (1979) and to Green and Godfrey 
(1983), this replacement can lead to common emitter current gain 3 * 100 in 
the traditional metal-contact bipolar transistor. 

This ten-fold increase in 3 corresponds directly to a ten-fold decrease 
of the recombination current in the n + emitter. Thus arises the implications 
for solar cells. The ten-fold decrease can only come from the presence of an 
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defective surface recombination velocity S at the n + Si/n + polySi interface. 
Recent work by NEUGROSCHEL (1984), outlined below, helped to some degree by 
the principal investigator, resulted in the conclusion indicated below 
concerning S. 

The physical mechanisms responsible for the decrease in S are complicated 

by a lack of knowledge of the minority-carrier diffusivity in both the n + 

monoSi and polySi , and by a poor present understanding of the interfacial 

layer. The parameters of this layer are highly sensitive to fabrication 

conditions and to surface treatment before the CVD process of polysilicon 
deposition. Various investigators have explored the chemistry of the 

interface, where peaks in P or As can occur, where deep-’evel impurities and 
interstitial oxygen may occur, and a thin insulating interfacial layer can be 
created by thermal oxidation or chemical treatments before CVD deposition. 

But one thing is clear. The polysilicon/silicon heterojunction acts as 
an effective surface passivant, reducing recombination losses. Various 

investigations suggest strongly that the ability to passivate persists in the 
presence of sunlight. 

Several possibilities exist for exploiting this passivant for Si solar 

i 

cells. Details about these possibilities appear in the recent proposal to JPL 

) 

by the author. 

We now sketch the measurement of S mentioned above. Consider an 

n + polySi/n + Si/p/n bipolar transistor. The n + Si emitter has a thickness of 
2xl0" s cm and a position independent doping concentration ~ 10 ^ As atoms/cm^. 
The n + polySi has a doping concentration of 5x10^ As atoms/cm'*, determined by 
SIMS, and a thickness exceeding 300 A. On control samples the n + polySi is 

l 

omitted; an ohmic contact is, made to the n + Si . In tne forward-active mode, 
the transistor shows a pre-exponential collector current of 10” 10 A/cm 2 (± a 
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factor of 2 ) on six samples and six controls, deriving from minority-electron 
injection into a base region having an impurity concentration ~ 10 ^ poron 
atoms/cm^. The base current and its pre-exponential factor, Jgj is also 
measured. 

Because Jgj derives almost entirely from holes injected into the n + 
emitter, we solve the minority hole continuity equation for this region. This 
yields Jg^ as a function of hyperbolic trigonometric functions having 
arguments involving the following characteristic velocities: 

(a) D E /W E = diffusion velocity, 

(b) W e /t e = volume recombination velocity, 

(c) S = surface recombination velocity at the n + polySi/n + Si interface. 

A standard expansion then yields 

J B0 * (en i /N eff )CCD E /W E )1 + + (V’E )] • (31) 

Here is the effective concentration in the n + Si , which takes into account 
Fermi statistics and energy-gap narrowing according to the experimental study 
of Neugroschel , Lindholm and Pao [18] and consistent with the theoretical 
model of Landsberg, Neugroschel, Lindholm and Sah [19]. The result is subject 
to the assumption of a parabolic quantum density of states as suggested by the 
experimental considerations of Keyes [20]-[21] and of Shibib and Lindholm 
[22]. Here also x t is the transit time for minority holes in the emitter 

t t /x E = (W^/2 D e x e ) + (W e /St e ) . (32) 

I 

I 

If we make the approximation, = t(Auger), then 
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( 33 ) 


1 + (x t /x E ) - 1 r W E /t(Aager)[(W £ /20 E ) ♦ (1/S) 

i -5 

* 1 * 40n ! {10 + 1/S} 

1 

■ 1. 

if S >> 400 cm/s, a conclusion that still holds if < x(Auger). 
Thus, 

°B0 (poly) = f l (U E /W E* S) 


J B0 (O > ‘ f 2< D E /M E> 

where J m ,(Q) denotes the pre-exponential current density when 
dU 

contact replaces the polysilicon. From the foregoing, 

S - {D E /W E ){CJ B0 CQ)/J B {poly)3-l) _1 • 

From the measured performance 
S * 5000 cm/s. 

Thus the polySi/N + Si interface provides an ohmic contact for 
electrons but a recombination velocity about three orders of magni 
than an ohmic contact for the minority holes. This supplies ample 
to investigate such heterojunction contacts for solar cells. 


(34) 


V35) 

an ohmic 
(36) 


majori ty 
tude less 
incentive 
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10. RECOMMENDATIONS ANO CONCLUSIONS 

I 

The short-circuit current decay method is apparently the only reliable 

i 

transient electrical method jfor determining x and S for a wide variety of Si 
solar cells. For very thin or very thick devices, this method must be 

supplemented. Otherwise only x or only S can be determined accurately, not 

both. The small -signal admittance method appears to fill this need. Our 

f 

i 

continuing work will explore this issue more fully. 

As recommendations, we offer the following: 

(a) Other laboratories, inlcuding JPL, would benefit by setting up the 
SCCO method. 

(b) The open-circuit voltage decay needs development for use in 

connection with manufacturing lines (see quarterly 2). This 

requires implementation of expressions for junction capacitance 

under forward voltage. These have been developed under this 

contract support, but are not yet written fully. 

(c) Measurement of the lifetime and the front surface recombination 
velocity of the quasineutral emitter is essential. An extension of 
the admittance methods used here may enable such measurements. 

(d) The polysilicon/silicon heterojunction has demonstrated properties 
desirable enough to warrant investigation about its use in solar 
cells. 
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APPENDIX A 


Although there are several ways to treat the sudden application of a 
short circuit replacing forward bias V, perhaps the simplest Is to think of 

i 

voltage -V being applied In series with V at t=0. This treatment emphasizes 

| 

the change in voltage that starts the ensuing transient. See Fig. A-l(a), 
bel ow. 

Thus at t 5 0, this change in voltage raises the right ohmic contact by 

i 

magnitude eV relative to the -left ohmic contact because the ohmic contacts are 
In equilibrium with the adjoining semiconductor in the sense that the distance 
between the quasi-Fermi level of majority carriers and the majority-carrier 
band edge remains the same as in equilibrium. They are in non-equilibrium in 
the sense that charge carriers can pass through the contacts. At t = 0 + , some 
arbitrarily small time after the application of the short circuit, the change 
in applied voltage has caused electrons to exit the n-type material adjacent 
to the contact, leaving behind unbared donor atoms and the positive charge 
shown in Fig. A-l(b). Similarly holes exit the p-type material (electrons 

enter the valence band from the metal), giving rise to the negative charge 
shown in Fig. A. 1(b). A near delta function of current i(t), flowing in the 
direction shown in Fig. A. 1(a), establishes this charge configuration at 
t = 0 + . Note that 1 ( t ) during the entire transient for t > 0 flows in a 

direction opposite to that occuring for negative time because the transient 
results in removing the electrons and hole present under forward V. 

Having established the existence of this negative charge, we now consider 
what happens subsequently. Here enters a result developed earlier from 
operating on Maxwoll's equation by the divergence operator: 

U = div curl H =di v j N + div j p + div[3(eE)/3t] . (A.l) 


5i 
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Figure A-l 

(a) For t < 0 , switch Si Is closed, So is open; conversely for t > 0; the 
junction space-charge region Is defined by x < x < x N . (b) charge density 
at t * 0 + . (c) charge density for t of the^order of a dielectric relaxation 
time, (d) the total current Is x-lndependent but is essentially majority- 
carrier convection current In the two quasi-neutral regions and is 
displacement current in the space-charge region for t of the order of a 
dielectric relaxation time, (e) electrons and holes drift out of the space- 
charge region in a transit time, (f) the resulting excess hole density In the 
space-charge region after a transit time has lapsed. 
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.From this result, two consequences emerge: (1) the charges In 

Fig. A-l(b) redistribute to the positions shown In Fig. A-l(c) within 

I 

the order of a dielectric Relaxation time x » c/o; and (2) the total 

! 

current Is solenoldal, that Is, Its divergence Is zero, where here the 

I 

total current Includes the displacement current. 

The consequence of (2) Is illustrated In Fig. A-l(d) for a 
particular time of order of t. Notice the large time-rate of change of 

I 

electric field E within the function space-charge region, x p < x < x N . 
Here we have employed a one-dimensional model so that the operator div 
becomes the operator 3/3x. Thus we see that the electric field in the 

space charge region grows rapidly so that within t of the order of t the 

barrier height has returned to its near equilibrium value and the 

electric field is several times larger than it was In negative time. 
But in negative time, the drift and diffusion tendencies of the junction 
space-charge region were perturbed only by perhaps one part in 10^ in 
the forward voltage steady state; that is, the space-charge region was 
in quasi-equilibrium. For t of the order r * t, the drift tendency now 
overwhelms the diffusion tendency, and holes and electrons drift out of 
the space charge region in a transit time x* determined by 

x^xp/velocity where the velocity approaches the scatter limited 
velocity because of the high field (Fig. A-l(e)). Hor typical devices, 
x' will be of the order of 10" 11 s. After this time has passed, the 
hole and electron concentrations will have returned nearly to their 
equilibrium values. 

Because 10"^ s is a time not observed by typical measurement 
equipment, we think of the initial condition established by shorting the 
terminals suddenly as that of quasi -neutral regions still storing 
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approximately the same excess charge as was present in the steady state 

I 

of negative time, of a spade-charge region at the equilibrium barrier 

I 

heiyht and a„void of excess holes and electrons, and as an excess 

i 

minority carrier density in the quasi-neutral regions that drops sharply 
to zero at the space-charge region edges (Fig. A-l(f)). 

I 

I 

APPENDIX B 

Determination of the Eigenvalues for SCCD and OCVD 
In this report, we have two Eigenvalue equations, Eqs. (9) and 
(12), that determine s^ of each mode for OCVD and SCCD. These are 

coth(X qN3 (1 ♦ V p ) 1 / 2 /L p ) + V 1 + s i Weff * U l9) 

and 

1 + (D (1 + s.x ) 1/2 /L S „)coth(X nMll (l + s. x ) 1/2 /L ) = 0 . (12) 

' p v l p' p eff' l)N8' l p' p' ' 

In (9) and (12), Eigenvalues exist only if 


1 + s i T p * 0 (° r s i < " 1/t p or x i > x p) 

where x. = -1/s . . - 

l l 

Granting (B-l), we have 


(1 


+ Vp 


) 1 / 2 =j(-l - Vp ) 1/2 


(B-l) 
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where (-I-s^t^) > 0. Replacing (1 + s^t in (9) and (11) with 
j ( - 1 - s i T p )1/2 yields 

cot ( X QNB K i /L p } - {D P K 1 /L P ) = ° (B ' 2) 


and 


1 + < D p K i /L p S eff )cOt(X 0NB K i /L p> = 0 


(B-3) 


where 


Equations (9) and (B-2) are identical and so are (11) and (B-3) 
under the condition of (B-l). (B-2) and (B-3) imply an infinite number 
of Eigenvalues as shown in Figs. B-l and B-2. 

For SSCD the vanishing determinant of the inverse of matrix A of 
Eq. (6) provides an alternate method for determining the Eigenvalues s^, 
but no such systematic method exists for OCVD. More details regarding 
this are planned for future publications. 



APPENDIX C 


To obtain the bbse diffusion length L and effective surface 

I 

recombination velocity S e ff at the back contact we have to solve the 

following two equations for the first-mode decay [2]: 

! 

i 

1 + L D pVVeff )c 0 t( V K l /L p) ' 0 < C1 > 

I 

I 

„ 01 qytWfrl 1 C0t[KiW QNB /Lj- [OpVL S eff ) 

1 P ( T p/ 2K i) + ( W QNB/ 2S eff ^ cosec ^ K i W QNB^ L p^ 

(C2) 

Here Kj = (-1 - SjTp)^^, = -I/tq-j, Wg^g is the width of the quasi- 
neutral base, and P(0, 0”) = (n?/N D p)[exp(qV/kT) - 1] where V is the 
steady forward voltage applied for t < 0. 

APPENDIX D 

The small-signal quasi-neutral base capacitance and conductance are 
given by [4]: 


,LF 

'QNB 


KD n 

*~ n 
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W p D n L nVeff 

L n ' °n 


S eff L n 


,W_ D 


(sinh^)(-jA:oth^- + S^)* 
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where K 


+ T 


D W 

C + S eff cotl> U £ 


n 


n ~TT 
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lj-— cot + S 


eff 


(Dl) 


>LF 

J QNB 


D W 

KD_ : "l” * S eff coth L 


n | n 
:U 


W 


" -Hcoth^ + S 
n L n 


eff 


(D2) 


= Aq(q/kT)(n?/N AA )exp[(qV/kT) - 1], 
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